Macroscopic traffic models from microscopic car-following models 
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We present a method to derive macroscopic fluid-dynamic models from microscopic car-following 
models via a coarse-graining procedure. The method is first demonstrated for the optimal veloc- 
ity model. The derived macroscopic model consists of a conservation equation and a momentum 
equation, and the latter contains a relaxation term, an anticipation term, and a diffusion term. Prop- 
erties of the resulting macroscopic model are compared with those of the optimal velocity model 
through numerical simulations, and reasonable agreement is found although there are deviations in 
the quantitative level. The derivation is also extended to general car-following models. 



I. INTRODUCTION 

For more 50 years, traffic flow has been a subject of 
intense research effort [0. While earlier studies were 
mostly conducted by traffic engineers, in the last decade 
the traffic flow problem has received great attention from 
the physics community as well, largely due to the sem- 
inal works dHl] m ^ ne early 90s, which demonstrated 
that traffic flow can be regarded as a driven nonequilib- 
rium system. There are empirical indications of multiple 
dynamic phases in the traffic flow and dynamic phase 
transitions . Several theoretical explanations 1 10 lf| 
for the empirical results were suggested. Also physical 
phenomena such as self-organized criticality and hystere- 
sis [jl6| were revealed. 

Numerous traffic models have been investigated (see 
Refs. Jl~T|— |l9|| for recent reviews) in relation to empirical 
data, and considerable progress has been achieved toward 
an understanding of various traffic phenomena observed 
empirically. Depending on the mathematical formulation 
used, traffic models may be categorized into one of the 
following types: car-following models, particle-hopping 
models, coupled-map lattice models, gas-kinetic models, 
and fluid-dynamic models. The first three types use a 
microscopic approach while the last type uses a macro- 
scopic one. The approach used in the gas-kinetic models 
is intermediate and may be called mesoscopic. 

Recently it was suggested [ 20 2lJ that different types of 
traffic models may belong to the same "universality" class 
in the sense that they share qualitatively similar prop- 
erties. More recently, a nonlocal fluid-dynamic model 
was derived from a gas-kinetic model p2| . These reports 
motivate further studies on mutual relationship between 
different types of traffic models. 

In this paper, we address the relationship between 
microscopic car-following models and macroscopic fluid- 
dynamic models. Specifically we use a coarse-graining 
procedure (Sec. |l|) to derive a macroscopic model 
(Sec. Ill) from the microscopic optimal velocity model, 
a particular case of the car-following-type model. The 
resulting macroscopic model consists of the continuity 
equation [Eq. (j^)] and a momentum equation [Eq. (Pq)] . 
The momentum equation contains a relaxation term, a 
density gradient term, and a diffusion term, similar to the 



fluid-dynamic model proposed in Ref. It is shown 
that both the density gradient term and the diffusion 
term arise from a directed influence due to the break- 
down of the balanced action-reaction. This is contrary 
to heuristic derivations |J] , in which the density gradient 
term is attributed to the velocity variance. It also pro- 
vides an origin of the diffusion term assumed in many 
fluid-dynamic models. In Sec. [V, the derivation is ex- 
tended to general car-following-type models. In Sec. [v|, 
the macroscopic model derived from the microscopic op- 
timal velocity model is examined numerically in compar- 
ison with the optimal velocity model. Section ^ con- 
cludes the paper. Some technical details are presented in 
Appendixes |a|, [b], and |c]. 

We remark that a different scheme to construct macro- 
scopic models from microscopic car-following models was 
proposed recently p^ |. The macroscopic fields p and 
v are defined via an interpolation procedure instead of 
a coarse-graining procedure. The resulting momentum 
equation is nonlocal, while our momentum equation is lo- 
cal. Also p and v defined in such a way do not strictly sat- 
isfy continuity equation (||) , while the continuity equation 
is an exact identity in the coarse-graining-based scheme. 



II. GENERAL FORMULATION 

In order to derive macroscopic traffic equations from 
microscopic ones, we first introduce two microscopic field 
variables, density field p{x,t) and flux field q(x,t), 



p(x,t) = 225{y n (t) - x) , 

n 

q(x,t) = 22y n (t)5 (y n (t) - 



(1) 



where y n (t) is the coordinate of the nth vehicle at time 
t with yi < y 2 < ■ ■ ■ < y n -i < y n < y n +i < ■ ■ ■ When 
traffic dynamics does not depend on third or higher or- 
der time derivatives of y n (t), these two fields specify the 
status of traffic flow completely. 

A natural way to obtain macroscopic description is 
to coarse grain these fields. We introduce a coarse 
graining envelope function <fi(x, t) which is non-negative 
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valued, peaked at (x, t) = (0,0), and normalized as 
J dxdt<t>[x, t) = 1. The coarse grained density p(x, t) and 
flux q(x, t) can be denned as 



p(x, t) = J dx'dt' <f>(x -x',t- t')p(x',t') , 

q(x, t) = / dx'dt' 4>{x -x',t- t')q(x', t') . (2) 



These two coarse grained fields specify the macroscopic 
status of traffic flow. 

Next we derive equations that govern the time evolu- 
tion of p(x,t) and q(x,t). For the evolution of p(x,t), 
one finds 



d . . d . . 
dt + dx q ^ X ' = 



(3) 



which describes the local conservation of vehicles in the 
coarse-grained description. This equation can be verified 
from Eq. (||) using integration by parts and change of 
variables. 

Derivation of the dynamic equation for q{x 1 t) is less 
straightforward. After some algebra, one obtains 



d d r 

— q[x, t) = p(x, t) (y n (t')) (Xtt) - [p{x, t) {yl{t')) [x t) 

(4) 

where the bracketed average of a quantity O n {x' , i') is 
defined as follows; 

(O n (x',t')) {xt) = — ^ J dx'dt 1 (/>(x-x',t-t') 

xJ20n(x',t')5(y n (t')-x'). (5) 

n 

Note that x', t', and n inside the brackets are dummy 
variables, while the label (x, t) in the subscript of the 
bracket notation represents a spatiotemporal position 
where the average is evaluated. This label will be omit- 
ted in the rest of the paper when its omission does not 
cause confusion. 

Here it is useful to introduce another macroscopic field 
v{x,t), 



v(x,t) = (y n (x',t')) =q(x,t)/p(x,t), 



(6) 



which represents some kind of macroscopic velocity, 
whose precise meaning depends on <f>(x,t). Two partic- 
ular coarse graining schemes are good for illustration: 
spatial coarse graining cf>(x, t) — 5(t)Q(X/2 — \x\)/X and 
temporal coarse graining (j)(x,t) — 5(x)@(T/2 — \t\)/T, 
where Q(x) is the step function which is one for x > 
and zero for x < 0. For the spatial coarse graining, v(x, t) 
becomes 



v{x, t) 



where the primed summation runs over the vehicles in 
the range (x — X/2,x + X/2) at time t. The denominator 
is equal to the total number of vehicles within the range 
and thus v(x,t) represents the arithmetic mean velocity. 
For the temporal coarse graining, on the other hand, it 
can be verified that 



1 = J2' n [yn(tn(x))}- 1 

u O>*) Eli 



where the primed summation now runs over the ve- 
hicles that reach the point x within the time interval 
(t — T/2,t + T/2), and t n (x) represents the time at which 
the nth vehicle reaches the position x. Here y n (t) > 
is assumed. Thus v(x, t) represents the harmonic mean 
velocity measured at local detectors. 

It is straightforward to rewrite Eq. (||) in terms of p 
and v instead of p and q. Also expressing Eq. (Q) in terms 
of p and v, one obtains 



' dv dv \ d 

Pll^+V— )=p(y n ( t '))-—(p9), 



dt dx J 



dx 



(7) 



where 



9(x,t) = (y 2 n (t'))-v 2 (x,t) 



measures the degree of microscopic velocity variation. 
Note that the left-hand side of Eq. (0) corresponds to 
the total derivative Dv/Dt = dv/dt + vdv/dx. Thus the 
two terms on the right-hand side can be interpreted as 
macroscopic force densities. The first term corresponds 
to the coarse-grained average of microscopic "forces" that 
act on each vehicle. The second term, on the other hand, 
arises from the coarse graining itself. In equilibrium sys- 
tems, 6 is proportional to the local temperature, and the 
second term represents the force due to thermal gradient. 

The remaining job is to express the force terms in 
terms of p and v. However, it is well known that a rig- 
orous treatment of the force terms generates an infinite 
sequence of dynamic equations. Thus we instead develop 
approximations of the force terms in Sec. Ill, so that 
Eqs. (||) and (?]) form a closed set of equations. This 
scheme is partly motivated by the absence of empirical 
indication that the dynamics of the forces is important. 

A procedure to derive a macroscopic model is illus- 
trated for the optimal velocity model in Sec. Ill and for 
general car-following models in Sec. [V. In both sec- 



tions, traffic states are assumed to be almost homoge- 
neous. In this linear regime, products of differentiated 
quantities such as Ylm=i(^ lm< ^rn/dx lm ) become progres- 
sively smaller as M increases, where l m are integers and 
O m are arbitrary functions of p and v. Therefore, it is 
sufficient to retain terms with M = or 1 only, which 
simplifies the construction of a macroscopic description 
considerably. In this sense, terms with M = or 1 can 
be called linearly relevant terms, and terms with M > 2 
linearly irrelevant terms. Properties in the linear regime 
such as the dispersion relation for small amplitude waves 
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depend on linearly relevant terms only. Effects of the 
linearly irrelevant terms with M = 2 are discussed in 
Appendix |A|. 



III. OPTIMAL VELOCITY MODEL 



We first study the optimal velocity model 



yn(t) = \[V op (Ay n {t))-y n (t)}, 



(8) 



where the constant A represents a driver's sensitivity and 
Ay„ = y n +i —y n is the coordinate difference between the 
vehicle n and its preceding vehicle n + 1. V op (Ay) is the 
optimal velocity to which drivers want to adjust their 
speed. An example is lop (Ay) = tanh(Ay — 2) + tanh2 
used by Bando et al. j24j. Here we will assume neither 
a particular functional form for V" op (Ay) nor a particular 
value for A p5[ . 

The coarse graining of Eq. (JsJ) leads to 

(y n ) = X[(V op (Ay n ))-v]. (9) 

The expansion of (Vo P (Ay n )) with respect to (Ay n ) gives 

(V op (Ay„)} = £ — ((Ay„)) ((Ay n - (Ay„))™} 

oo 

= V op «Ay n » + Im, (10) 

m— 2 

where I m is the term that is proportional to 
((Ay„ — (Ay„)) m ). Here I\ is absent since 

((Ay„ — (Ay„))) — 0. Note that the leading correction I 2 
compensates for the difference (Vo P (Ay„)} — V op ((Ay n )), 
which is positive (negative) when V op is a convex (con- 
cave) function. In the linear regime, however, all cor- 
rections I m (m > 2) can be ignored. Moreover it can 
be shown that the second term on the right-hand side 
of Eq. (ph is also negligible in the linear regime (see Ap- 
pendix B). Therefore, the derivation of a macroscopic 
description in the linear regime is reduced to developing 
a proper approximation of (Ay n ). 

A. Directed influence 

A naive approximation of (Ay„)( x ( ) is p^ 1 {x,t). How- 
ever, this seemingly reasonable approximation has a se- 
rious problem. For illustration, it is useful to introduce 
an unphysical model by replacing Ay n {t) in Eq. (^) with 
Ay„-i(t), so that each vehicle responds to the vehicle 
behind it rather than the vehicle ahead of it. This un- 
physical model, which differs from the physical one only 
by the directionality of the influence, has qualitatively 
different properties. Thus proper macroscopic descrip- 
tions should contain information about the directionality, 



while a naive approximation fails to capture this infor- 
mation. 

To take the directionality into account, an intuitive 
prescription was proposed |26| without a rigorous justifi- 
cation, 



< A JM>(w) ^p- X {x + l/2p{x,t),t), 



(11) 



which amounts to evaluating the density at the midpoint 
between two vehicles n and n + 1. For the above un- 
physical model, this prescription results in an expression 
which is similar to Eq. (fill) but has a negative sign in 
front of 1/2. Thus this prescription contains information 
about the directionality. 

In the linear regime, we find that a controlled approx- 
imation of (Ay„) can be obtained in a rigorous way (see 
Appendix £j ■ The result is 



(Ay„) 



1 dp- 1 
2p dx 



(12) 



where £ represents the sum of all terms with second 
or higher order derivatives. Note that Eq. ( |l2| ) agrees 
with the Taylor expansion of the heuristic approximation 
[Eq. (pH)]) up to the first order derivative correction to 
p^ 1 . The deviation occurs in the second order derivative. 
While the second order derivative in the Taylor expan- 
sion of Eq. ( ]i"l| ) comes with the coefficient 1/8, a rigorous 
calculation leads to the coefficient 1/6 (see Appendix^): 



E = 



1 d 2 p- x 



O 



dV 1 



(13) 



6p 2 dx 2 \ dx 3 

Thus the leading term in Eq. (|l^) can be expanded as 

1 dp- 1 



V op ({Ay n )) = V op (p~ 1 ) + V'(p- 1 ) 



2p dx 



(14) 



where Ei r denotes the sum of linearly irrelevant terms. 
By combining Eqs. (|?|), (||), (^0|), and (|T4|), one obtains 



dv 
di 



dv r 

v— = A [V OJ 



dx 



2p 



K 



dp- 



dx 



(15) 



Note that the second term proportional to the density 
gradient arises from the directed influence, while conven- 
tional derivations of fluid-dynamic models Q] attribute 
the density gradient term to the velocity variance term 
in Eq. (Q) . We will call the second term the anticipation 
term. The first term is often called the relaxation term. 

It is interesting to compare the dispersion relations of 
microscopic and macroscopic models. In a microscopic 
description, small perturbations with respect to the ho- 
mogeneous state can be written as 



y n (t) =v h t + p h 1 n + 6yexp(iKn + 'yt), 



(16) 
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where Vh — V op (p h 1 ) . By linearizing Eq. 
the dispersion relation 



7± 



-1 ± A/1 



A 1 



1) 



one obtains 



(17) 



On the other hand, small perturbations in the macro- 
scopic description can be written as 



p(x,t) = ph + 8pex^>(ikx + u>t), 
v(x,t) = Vh + Sv exp(ikx + tot), 



(18) 



where kpj^ 1 is the macroscopic counterpart of n since 
both represent the phase difference between two succes- 
sive vehicles, and tu + ikvh is the macroscopic counterpart 
of 7. To see the origin of the additional term ikvh, note 
that 7 is the frequency measured in the moving reference 
frame with the velocity Vh , while to is the frequency mea- 
sured in the stationary frame. By linearizing Eqs. (|^) 
and (|l5|), one finds 



where 



ikvh = — 



-l±\ 1 



op 



Mkpl 1 ) 



A(x) — ix + 



(ix) 2 



(19) 



(20) 



when the last term in Eq. (15) proportional to E is ig- 
nored and 



A(x) = ix + 



(ix) 2 (ix) 3 



(21) 



2 6 

when the leading contribution to E in Eq. (|l3| ) is in- 
cluded. Note that A(x) agrees with the Taylor expansion 
of the factor (e iK -l) in Eq. @. Thus it is clear that the 
macroscopic momentum equation ([l5|), combined with 
the continuity equation (|^) , gives a correct description of 
the long wavelength behavior of the microscopic model 
[Eq. (0)] in the linear regime. 



B. Effective diffusion 

Despite the excellent agreement of the long wavelength 
components, it is premature to accept Eq. (jl^) as a 
macroscopic momentum equation since naive treatments 
of £ introduce an artificial instability, which is absent in 
the microscopic model [Eq. (g)]. For demonstration, we 
examine the linear instability criteria. In the microscopic 
model, from Eq. (jlT]) one obtains that small fluctuations 
of the mode k become linearly unstable when 



K' P (Ph 



> 



x 



1 



(22) 



Note that the n = mode shows the strongest instability 
and at the critical density where the instability first sets 
in, only an infinite wavelength mode becomes unstable. 



In contrast, naive macroscopic models give different re- 
sults. When S is ignored completely, Eqs. ([l9]) and j20| ) 
result in a linear instability criterion V^p^ 1 ) > A/2 
for mode k. Note that this inequality does not contain 
k. Thus as soon as ph satisfies this inequality, fluctua- 
tions of all wavelengths become unstable simultaneously, 
different from the behavior in the microscopic descrip- 
tion. On the other hand, when the leading contribution 
to E in Eq. (13) is retained, Eqs. (p^ ) and ( f2~l| ) result in 
KpiPh 1 ) > W2)[l - (^Ph 1 ) 2 /6]"^Note that the right- 
hand side vanishes as kp^ 1 — » 00 and thus the homoge- 
neous state is always unstable with respect to fluctuations 
with small wavelengths. This artificial instability cannot 
be cured by merely using higher order approximations of 
E. For example, if we assume that the next order con- 
tribution to £ is (l/4\p 3 )(d 3 p^ 1 /dx 3 ), which generates 
the correct next order in A(x), one obtains the linear in- 
stability criterion V^ 1 ) > (A/2)[l - (kp- 1 ) 2 / 12]/ [1 - 
(kp^f/Q] 2 , which again shows an artificial instability 
for the short wavelength components. 

To find the origin of the failure, it is useful to analyze 
the microscopic dispersion relation [Eq. ([l7])] since the 
approximations of E are equivalent to truncating the se- 



ries e 



1 



; + (ik) 2 /2 + (i K ) 3 /3\ + (m) 4 /4! + 



at 



a certain order. It can be verified that when the series is 
truncated at a, finite order, highest order terms dominate 
the physics for large k and generate the artificial insta- 
bility for large n (» 1) modes, while such instabilities 
are absent when the series is summed up to the infinite 
order. Thus it is clear that truncation at a finite order is 
responsible for the artificial instability. 

In this subsection, we aim to develop an approxima- 
tion of E, which is compact but still captures important 
features of the exact E. A key observation is that modes 
with kp^ 1 3> 1 are unphysical since fluctuations on length 
scales shorter than the vehicle spacing are not defined in 
the original microscopic model. Motivated by this ob- 
servation, we transform the leading order term of E in 
Eq. ( p"3| ) in such a way that it preserves the same long 
wavelength behavior but suppresses fluctuations in short 
wavelength components with k 3> Ph ■ To implement this 
idea, one first notes that Eq. (||) relates small fluctuations 
of p and v as follows: 



x ik P h x 
op = — ov. 

uj + ikvh 



(23) 



One then exploits the correspondence between to + ikvh 
and 7, and between kp~^ 1 and k. From the result 7+ «; 
V 

op 



Ph' 



for small n, one obtains 



Sp~ 



■ 5v. 



In this derivation, the 7_ mode is ignored since it al- 
ways decays with time. Note that the resulting rela- 
tion amounts to a variational form of v = V op (p^ 1 ) that 
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can be regarded as the zeroth order approximation when 
kp^ 1 C 1. Its first or higher order corrections will be 
ignored since they introduce third or higher order deriva- 
tives to the new approximation of £ [Eq. (E4|)] . This way, 
we construct an approximation 



i d 2 P - 1 



6p 2 dx 2 6p 2 dx 
The momentum equation becomes 



dv 
di 



dv 



AVJ P dp 



A d 2 v 



2p3 dx 6p 2 dx 2 



(25) 



Note that our approximation of £ results in a diffusion 
term, which tends to suppress short wavelength fluc- 
tuations. Indeed, the linear instability criterion from 
Eqs. |) and (H) becomes V> p (p^ 1 ) > A(l + k 2 /6p 2 ) 2 /2, 
which confirms the suppression of modes with k 3> Ph- In 
addition, it can be verified that the macroscopic and mi- 
croscopic dispersion relations agree up to order k 3 . Thus 
we conclude that Eq. (|2^) is a satisfactory macroscopic 
momentum equation in the linear regime. 

Finally, we remark for completeness that Eq. ( ^5| ) can- 
not be used to study backward time evolution. This re- 
striction arises from the neglect of the 7_ mode, whose 
magnitude does grow in the backward time evolution. 



IV. GENERAL CAR-FOLLOWING MODELS 



In this section, we extend the derivation in Sec. Ill 
to general car-following models. When third or higher 
order time derivatives do not appear in microscopic traf- 
fic equations, a general car-following equation with the 
Galilean invariance can be written as 



Vn = A op (Ay n ,Ay n ,y n ). 
Coarse graining leads to 

dv dv . , . . . , . . . . 
-Q t + « A op {(Ay n ), (Ay n ),v) , 



(26) 



(27) 



where (Ay n ) can be approximated by Eqs. fll2| ) and (13), 
and 



1 dv 1 d 2 v 

See Appendix for a derivation of Eq. (| 
expand A op (- as 

4>p(---) ^ A op {p-\0,v) 



(28) 
We further 



+ A. 
+ A op , 2 



f 1 dp- 1 1 d 2 p- 1 
''''' \2p dx 6p 2 dx 2 
ldv 1 d 2 v \ 
p dx 2p 2 dx 2 J ' 



(29) 



where A opA = 9 2s A op (zi, z 2 , z 3 )|( 2l , 22 , 23 )=(p- 1 ,o,i))- In 
real traffic systems, A)p,i and A oPt 2 are expected to be 
positive while A op ^ is expected to be negative. Cross- 
terms proportional to A op ^A op ^2 are ignored since they 
are linearly irrelevant. The macroscopic momentum 
equation then becomes 



dv dv , _, s 



A^dp- 1 A opA d 2 p- 1 



2p dx 



6p 2 



+ ■ 



lop, 2 



dv 



A> P ,2 d 2 v 



dx 2p 2 dx 2 ' 



dx 2 
(30) 



Note that the dependence of A op on Ay n gives rise to an 
explicit diffusion term. 

Despite the explicit diffusion term, the artificial in- 
stability at short wavelength components may still arise 
when A oPj i is sufficiently large since the term propor- 
tional to d 2 p~ 1 /dx 2 tends to generate the artificial in- 
stability, as demonstrated in Sec. III. Thus we follow the 
procedure in Sec. Ill B to obtain 



d 2 p~ 



dx 2 



A oPi3 d 2 v 
A oP: i dx 2 



(31) 



which is a generalization of Eq. ( |24| ) . The resulting mo- 
mentum equation is 

dv dv , _ t v A op ,i dp' 1 
-+ v 7iZ= A °p(p >M + - 



dt dx upyr ' ' > 2p dx 

. A op .2 dv 3A op , 2 - ^o P ,3 9 2 v 



p dx 



6p 2 



dx 2 



(32) 



Note that the factor 3A op ,2 — A op .^ in front of the dif- 
fusion term is manifestly positive. This equation is free 
from the artificial instability. 

To elucidate the relation with Eq. (p5|), it is useful to 
define an effective optimal velocity Vo PiC ff (p -1 ) in an im- 
plicit way as a solution of 



A op (p \0, V op ,eff) = 0. 



(33) 



When v4 op ,3 < for all v, the solution is unique and there 
is no ambiguity in V oPje g(p ). One also defines 



A e ff(p 1 ,v) 



Vop,Gff(P _1 ) - V 



(34) 



which is positive for all p and v if ^4 p,3 < always. 
Thus the first term in Eq. (|3|) can be interpreted as a 
generalized relaxation term: 



A op (p 1 ,0,u) = A e ff(/0 1 ,v)[V ap , e s(fi l ) 



(35) 



In certain situations, the third term in Eq. ( |32|) can 
be transformed into a familiar form. One applies the 
procedure in Sec. |l| B to the term, and uses the re- 
lation 7+ w — (A oPi i/Ao Pj a)iK(l + f3in), where (3 = 
1/2 - A oPy2 /A op ^ - A oP: i/Al p 3 . Thus we obtain 
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dv 
dx 



lop,! 



dp 



A> P ,3 dx 



f3 2 v 
p dx 2 



(36) 



where third or higher order derivatives are neglected. On 
the other hand, the second order derivative should be 
kept since it renormalizes the diffusion term. The macro- 
scopic equation of motion then becomes 



dv dv . rT , 
_ +u _ = A cff [K P „ 



vA opl dp pA oP:3 d 2 v 



2p 3 dx 



6p 2 



dx 2 ' 
(37) 




where v = 1 - 2A op , 2 /A op:i and p = 1 - 3A p i2 /A>p,3 - 
6/3A op . 2 /A op . 3 = 1 - 6(A p )2 / J 4op,3)(l - A p l2 /A>p,3 - 
i ). Note that three force density terms in 
are in one-to-one correspondence with those in 
Moreover the corresponding terms in the two 
equations usually have the same sign since v is positive 
and ^4 p,3 is negative. However when (3 in Eq. ( |36| ) is a 
sufficiently large negative number, p in the diffusion term 
in Eq. ( p7| ) becomes negative, and an artificial instability 
at short wavelength components arises. Thus Eq. ( |37| ) 
can be used only when p is positive while Eq. (|3^) can 
be used in general situations. 



V. MICRO VS MACRO 

In this section, we compare numerically the properties 
of the microscopic optimal velocity model [Eq. (JSJ) ] and 
the macroscopic model [Eqs. (||) and (]25|)] derived from 
it. For definiteness, we use 



V op (Ay) = 



tanh 2 



^neutral 
-^width 



with v max = 33.6 m/s, £ n cutrai = 25.0 m, a; W idth = 
23.3 m, c bias = 0.913, and A = 2 sec -1 as in Ref. p7j. 
A system size L = 2.33 km is simulated with N vehi- 
cles (ph = N/L), and the following microscopic initial 
conditions are used: 



np^ 1 + Asm(&nnp h /L), 1 < n < N/3, 



Vn{0) 

y n (Q) =np h l , 

y n (0) = V op (Ay n (0)), 



N/3 <n< 2N/3, 
for all n. 

(38) 



The corresponding macroscopic initial condition is pre- 
pared by coarse graining the microscopic initial condition 
[see Eqs. (|]) and (jfy] with the spatial coarse graining 
function 4>(x,t) = (27ro- 2 ) _1 / 2 exp(— x 2 /2a 2 )5(t), where 
we choose a = 46.4 m. The periodic boundary condition 
is imposed for both the microscopic and macroscopic sys- 
tems. 

We first verify that the density range p c \ < p < p c2 , 
in which the homogeneous traffic state becomes unsta- 
ble with respect to infinitesimal perturbations, is essen- 
tially identical for the microscopic and macroscopic mod- 



els. This implies that, in the linear regime, the macro- 
scopic model describes the long wavelength behavior of 
the microscopic model very accurately. 

To quantify the accuracy of the macroscopic model, 
we introduce the space-averaged relative deviation d v it) , 
which is defined by 



d v (t) 



\/ ([^macro(^i t) ^micro^) 0]^)spa 
(^micro(^> 0) space 



where (■ • -)spacc represents the spatial average. Here 
fmacro(£, t) is calculated from the macroscopic model, 
while v m icro(:£, t) is obtained by coarse graining the mi- 
croscopic configuration at the time t. 

When the initial perturbation from homogeneous flow 
is small, say A = 1.165 m, we find that d v (t) is negligible 
for all density outside the linearly unstable density range. 
A typical velocity profile is shown in Fig. |l} Note that the 
macroscopic profiles are almost indistinguishable from 
the microscopic ones. Even when N = 72 (131), which 
corresponds to a density slightly below (above) the lower 
(upper) critical density p c i[ c2 ) ~ 73 (130)/2.33 km (nu- 
merically obtained critical densities are nearly the same 
as analytic ones), d v (t) remains ~ 2x 10~ 4 during several 
hours of simulation time. 

The accuracy in the linearly unstable density range is 
also examined for A — 1.165 m and N — 73, which is the 
smallest N that demonstrates the linear instability. The 
microscopic simulation shows that the initially smooth 
profile becomes "rough" as short wavelength fluctuations 
develop. An almost identical roughening is found in the 
macroscopic simulation, and d v (t) is almost negligible ini- 
tially [Fig. |2|(a)]. However, the growth rate of the short 
wavelength fluctuations is faster in the microscopic sim- 
ulation compared to the macroscopic simulation. This 
difference is responsible for the rapid growth of d v (t) 
near t ~ 55 min. The growth of d v (t) occurs at an ear- 
lier time for the density with stronger linear instability. 
Both in microscopic and macroscopic simulations, after a 
sufficient time interval (< 120 min) all short wavelength 
fluctuations merge into a single large traffic jam, which 
moves backward at a constant speed without further evo- 
lution in its shape. Thus this jam corresponds to the final 
steady state. Figure |2|(b) compares the velocity profiles 
of the jams from the microscopic and macroscopic sim- 
ulations. The velocity of the jam propagation speed is 
different and the locations of the jams coincide periodi- 
cally in time, resulting in the periodic dips in Fig. ^|(a). 

Next we choose A = 74.56 m in Eq. d3cf), and exam- 
ine the performance of the macroscopic model for large 
perturbations. Figure ||(a) shows the initial density pro- 
file. After a sufficiently long time, the initial condition 
may evolve to a homogeneous state or to a congested 
state. The evolution to a congested state is realized 
for 65 < N p$ 156 when the microscopic model is used 
and for 66 < N < 147 when the macroscopic model is 
used. Thus the lower critical density is in good agree- 
ment while the upper critical density shows about 6% 
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deviation. The comparison with the linear critical densi- 
ties shows that both microscopic and macroscopic models 
exhibit metastability, which implies the hysteresis phe- 
nomena in the metastable density range. The phase di- 
agram in Fig. ^(b) summarizes the result. Note that the 
microscopic metastable regions are wider. 

We also investigate the dependence of the critical den- 
sity on A for fixed A = 74.56 m. It is convenient to in- 
troduce a dimensionless parameter A = (xwidth/^max^, 
which is about 1.387 for A = 2 sec -1 . Figure [|(a) shows 
the relative deviations of the macroscopic critical densi- 
ties with respect to the microscopic ones. For the lower 
critical density, the macroscopic result is in good agree- 
ment with the microscopic one for general A. For the 
upper critical density, on the other hand, the deviation 
of about 6% at A w 1.387 shrinks with the increase of 
A and good agreement is achieved near A = 2. Thus 
the difference between the microscopic and macroscopic 
metastable regions in Fig. |^(b) shrinks as A — > 2. 

The velocity — v g of a backward propagating traffic jam 
cluster (v g > 0) is also investigated. Since v g is almost 
independent of N, we fix N = 100 (p h w 42.9_km -1 ) 
for simplicity, and examine v g as a function of A. Fig- 
ure ^(b) (diamonds) shows the ratio between the micro- 
scopic value w™ lc and the macroscopic value w™ ac . Note 
that u™ IC /z;™ ac f» 1 when A is close to 2. This agreement 
is notable considering that the macroscopic model does 
not have any free parameter which can be varied to en- 
hance the agreement. The agreement, however, becomes 
less satisfactory as A becomes smaller. 

A crude understanding for the good agreement near 
A = 2 can be achieved via the linear analysis, although 
the given initial condition is not in the linear regime. For 
the general optimal velocity model, the linear instability 



(dv/dx)\ (dp 1 /dx) 2 , and (dp 1 / dx](dy I dx) (see Ap- 



develops when V' > 



A/(l + cos k); here we introduce 
This inequality sets an upper 



Kp = (Zwidth/Umax)Kp 
limit K c , above which the instability does not appear. 
Note that k c shrinks to zero as A/2 approaches max(V^ p ), 
which is 1 . Thus the characteristic length scale of the in- 
stability becomes longer as A — > 2. This may explain the 
excellent agreement near A = 2, since the macroscopic 
model becomes more precise as the characteristic length 
scale grows. 

From these comparisons, we conclude that the macro- 
scopic model [Eqs. (||) and (|||)] is quite accurate in the 
linear regime, and provides a reasonable description of 
fully developed jam clusters in the nonlinear regime, al- 
though there are deviations in the quantitative level. But 
when short length scale dynamics plays an important 
role, for example when the avalanchelike growth of many 
small clusters occurs, the macroscopic model is not sat- 
isfactory. 

To construct more accurate macroscopic models, one 
needs to take into account effects of various terms ignored 
in the macroscopic momentum equation derivation. As 
a first trial, we extend the derivation to the nonlinear 
regime by including effects of all terms proportional to 



pendix [AJ) . The resulting equation (|A8[ ) for the same 
optimal velocity model is examined. As expected, the 
linearly unstable density region is identical to that by 
Eq. (|2^). However, thejratio w™ lc /v™ ac deviates further 



i i' i-i y ' y 

from one [circles in Fig. |4j(b)]. Thus it appears that naive 
inclusion of linearly irrelevant terms does not improve the 
accuracy. 



VI. SUMMARY 

A local macroscopic fluid-dynamic model is derived 
from a microscopic car-following model, which establishes 
the link between the two types of traffic models. It is 
emphasized that the directed influence due to the break- 
down of the balanced action-reaction is an important 
ingredient. For the optimal velocity model, the corre- 
sponding macroscopic momentum equation consists of a 
relaxation term, an anticipation term (proportional to 
the density gradient), and a diffusion term. Thus it has 
a structure similar to the fluid-dynamic model in Ref . [Q . 
However, the density gradient term is found to arise from 
the directed influence rather than the velocity variance. 
It is demonstrated that the diffusion term also arises from 
the directed influence. The derivation provides an unam- 
biguous way to determine the coefficients of the anticipa- 
tion term and the diffusion term. The macroscopic model 
derived from the optimal velocity model is examined nu- 
merically, and its properties are found to be in reasonable 
agreement with those of the microscopic model although 
there are deviations in the quantitative level. 
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APPENDIX A: EFFECTS OF LINEARLY 
IRRELEVANT TERMS 



While the derivation in Sees. HI and |^ assumes a lin- 
ear regime, interesting traffic phenomena often occur in 
the nonlinear regime. In this appendix, we aim to develop 
a macroscopic momentum equation, which is applicable 
to nonlinear traffic phenomena when the characteristic 
length scale is sufficiently long. For traffic phenomena 
with a long characteristic length scale £, each derivative 
d/dx can be formally regarded as a small expansion pa- 
rameter since it effectively introduces the small factor 
l/£. Then we can take a perturbative approach: terms 
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without derivatives constitute the zeroth order contribu- 
tions, and terms with the first order derivative the first 
order contributions. Thus the relaxation and anticipa- 
tion terms are the zeroth and first order contributions, 
respectively. All zeroth and first order contributions are 
already included correctly in Eqs. (|2^) and (^). As for 
the second order contributions, however, only part of 
them are included since terms proportional to (dv/dx) 2 , 
(dp- 1 /dx) 2 , or (dp- 1 / dx)(dv / dx) are of the same order 
as the diffusion term. Below we demonstrate a procedure 
to obtain the missing second order contributions for the 
general microscopic model [Eq. (J26j)] . 

In the general expression (^), the last term propor- 
tional to d(p8)/dx is irrelevant for our discussion since 
it generates third or higher order contributions only (see 
Appendix |bJ). We expand the first term to obtain 

(y n (t')) ^A op ((Ay n ),(Ay n ),y n ) 
+ ^fl((Ay n -(Ay n )) 2 

A>p,22 



2 

A op ,33 



(Ay n - (Ay n )Y 



+ A opA2 ((Ay n - {Ay n ))(Ay n - (Ay n ))) 

+ ^o P: 23((Ay„ - (Ay n ))(y n - v)) 

+ A opA3 ((Ay n - (Ay n ))(y n - «)), (Al) 

which is a generalization of Eqs. (^) and (|l(i|). Here 

A op ,ij = d z A i A ov (zi, Z2, za)\( z , ■*».* H i=fo- 1 .Q.t>v I n 
Sees. [II and Evl, the last six terms in Eq. (Al) have 
been ignored. For a spatial coarse graining function 
4>(x,t) = 4>x{x)5{t), we find 



.4 



op, 11 



(Ay n - (Ay n )Y 



a 2 A 



op, 11 



A, 



op, 33 



(Vn - v)' 



a 2 A 



op, 33 



A op ^ 3 ((Ay n y n - (Ay n )v)) 
where a 2 = J dx'x' 2 (f>x(x'). 



(J 2 A op ,i3 



( dp- 1 
V dx 

Ox 
dv dp -1 



(A2) 



dx dx 

Note that these second or- 
der contributions depend on the coarse-graining func tion 
explicitly The other three nonlinear terms in Eq. (Al) 
give third or higher order contributions only (see Ap- 
pendix H) . 



The first term on the right-hand side of Eq. (Al) 
also generates the second order contributions. The sec- 
ond order expansion of its arguments results in (see Ap- 
pendix ^J) 



(Ay n 
{Ay n 



1 dv 
p dx 



1 dp' 1 1 d 2 p- x 
6p 2 dx 2 

v ( dp -1 



2p dx 
1 d 2 v 
2p 2 dx 2 



1 



dp- 1 
dx 



2 V dx 



(A3) 



Thus one finds 



A op ((Ay n ), (Ay n ), (y n )) 
A op (p-\0,v) 



A 



P .i 



I dp- 1 + 1 d 2 p- 1 



1 fdp- 



+ A oPi2 
. A opA i f dp 



2p dx ' 6p 2 dx 2 6p \ dx 

1 dv 1 d 2 v v ( dp- 1 ^ 2 
p dx 2p 2 dx 2 2 \ dx 

-1\2 , / Q .x2 



\ dx 



Ao P , 



p,22 



2p 2 



dv 
dx 



A 



op, 12 



(A4) 



dp 1 dv 



2p 2 dx dx 



Note that the second order contributions from the ex- 
pansion of A op ((Ay n ) , (Ay n ), (y n )) do not depend on the 
coarse-graining function. Next we apply the prescription 



d f dp 



dx V dx 



d_ 

dx 



A op ^ dv 
A op .i dx 



(A5) 



which is the extension of Eq. ( ^l|) to second order. The 
resulting macroscopic momentum equation is 



dv dv 
dt dx 
A op (p~\0,t;) 
A op ,i dp- 1 
2p dx 



A> P ,2 dv 3A oPj2 - A op ^ d 2 v 



2 dx 



6p 2 



+ 



A op ,i 


vA op , 2 


A op! n a 2 


6 P 


2 


8p 2 


A op .22 


^4op,33 


A)p,3Ajp,13 


2p 2 


6p 2 


6p 2 A op ,i 


A)p,12 


^.op,13 


A)p,3Ajp,ll 


2p 2 


6p 2 


6p 2 A opA 



iop,ll 



dx 2 

( dp- 1 
V dx 



o 2 A 



op, 33 



cr 2 ^4op,13 



(-) 

\dxj 
dp- 1 dv 
dx dx 
(A6) 



For the optimal velocity model [Eq. (g)], this reduces to 
dv dv _ x XV^dp- 1 A d 2 v 

'VL V'L a 2 V' r 



+ X [-6p- + 8p^ + ^~ 



dx 



A V" p dp 1 dv 
6p 2 V^'p dx dx 



(A7) 



From numerical simulations we find that the last two 
terms give rise to the artificial instability for short wave- 
length components despite the presence of the diffu- 
sion term. It turns out that the artificial instability 
can be cured by applying the prescription (dp- 1 /dx) ss 
(l/V/ ip )(dv / dx) . Thus the resulting momentum equation 
for the optimal velocity model reads 
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^Kp dp 1 A d 2 v one finds 



® V dv hi, / -In i 

- + V - = A[T4 p (p )-«] , 2p Qx , 6/?2fe2 



(A8) 



APPENDIX B: IRRELEVANT TERMS IN THE 
LINEAR REGIME 

In this appendix, we assume the spatial coarse graining 
4>{x,t) — 4>x(x)5{t) for dcfinitcncss. 

(i) 9 = ((y n — (y n )) 2 )' After some algebra, it can be 
written as follows, 



9(x,t) 



2p 2 



dx' dx <f>x(x — x)4>x(x — x ) 



When the characteristic length of the variations is much 
larger than the coarse-graining scale, m — n can be for- 
mally regarded as small numbers. To obtain the leading 
contribution, we may then use the formal approximation 



Vm(t) - y n (t) 
which leads to 



dv 
dx 



[x,t) 



[Vm{t) - y n (t)}, 



Note that the second factor on the right-hand side is 
proportional to the square of the spatial extension of the 
coarse-graining function. When there are many vehicles 
within the coarse-graining scale, 

{yl)-{y n fno\ 

where a 2 = J dx'x' 2 (f>x(x'). Thus we obtain 



(x,t) fa a 2 



dv 

dx 



(ii) ({Ay n — (Ay n )) 2 ): The procedure is very similar: 
((Ay n -(Ay n )) 2 ) (xt) 

= — r / dx'dx"(j)x(x — x')d>x(x — x") 
2p 2 J 

x [Ay m (i) - Ay n (i)] 2 S(y m (t) - x')S(y n (t) - x"). 



Using the formal approximation 



Ay m {t) - Ay n (t) 



dp 
dx 



■[vm{t) - y n (t)}, 



((Ay n -(Ay n )) 2 )*a 2 

(hi) ((Ay n - (Ay n )) 2 ): 

((Ay n -(Ay n )) 2 ) (xt) 
= J dx'dx"cj>x{x-x')cl)x{x-x") 
x [^Vm(t) - Ay n (t)} 2 S(y m (t) - x')8{y n {t) - x") 



Since {Ay n ) w (l/pVdv/dx) in the leading approxima- 
tion (see Appendix. |^) , we use the formal approximation 

Ay m {t) - Ay n (t) » A [y m (t) - !/»(*)], 



and obtain 

((A 2 /„-(Ay„)) 2 )«a 2 

(iv) <(Ay„ - (Ay n ))(Ay n - (Ay n )) 
((Ay n - (Ay n )){Ay n - (Ay n )) « a 

(v) ((Ay n -(Ay n ))(y„-w)) 



1 2 



d_ (ldv 
dx \p dx 



2 dp 1 d f 1 dv 



dx dx \p dx 



((Ay n -(Ay n ))(y n -v)) « a 2 



(vi) ((Ay n - (Ay n ))(y n -v)): 



d ( 1 dv \ dv 



dx \p dx J dx 



((Ay n - (Ay n ))(y n - v)) « a 



2 dp 1 dv 
dx dx 



APPENDIX C: MACROSCOPIC EXPRESSIONS 
FOR THE DIFFERENCES 

This appendix presents derivations of Eqs. (|l^), (|l2|), 
and @. 

(i) (Ay n ): One begins with the definition of (Ay n ): 



p(Ay n ) — J dx dt 1 4>{x — x' , t — t') 

x E^+i^') - *»(*')]*(»»(*') - *')• (CI) 

71 

The following identity is useful: 

J2 [Vn+i(t) - y n (t)} 8(y n (t) - x) = —y r(Xtt) (t), (C2) 
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where r(x, i) is the vehicle number right in front of x 
at time t. For example, when y m {t) < x < y m +i(t), 
r(x,t) = m + 1. Note that each side of the equation 
vanishes unless there is a vehicle at x, and that the x in- 
tegration of each side from y m (t) — e to y m (t) + e results 
in y m+ i(t) — y m (t) , which proves the identity. Using the 
identity, Eq. ( |C l| ) can be simplified to 

(Ay n ) = p- x + p- x 2-[Ax(x,t) + A 2 (x,t)], (C3) 



where 



A 1 (x,t) = J dx'dt'cj)(x - x',t -t') 

y r (x',t')(t') - y r ( X ' ,t')-i(t') 



l d P - 2 (x,t) p-Hx.t) dp-^t) 

A2(x ' t} "l2 dx =— ■ (C7) 



dx 



From Eqs. (E3) and (C7), one then finds 



1 d 2 p" 1 1 /dp" 1 
6p 2 dx 2 6p \ <9x 



(ii) (Ay n ): For derivation, it is convenient to relate 
(Ay„) to (Ay n ). Using integration by parts, one can 
verify 

p(Ay n ) = — [p{Ay n )] + — [p(y n Ay n )} . (C8) 



A 2 (x, t) = j dx'dt'4){x -x',t- t 1 ) 

~y r (x',t')(t') + y r (x',t')-i(t') 



(C4) 



To obtain Eq. (|C3|), the integration by parts is used. 
Below </>(x, t) is assumed to be even in x. In the ho- 
mogeneous state, y r (x',t')(t') - y r (x>,t>)-i(t') = p^ 1 and 
A\{x, t) = l/2p(x, t) since r(x',t') — 1 is the vehicle num- 
ber right behind the position x' at time t' . It can also be 
shown that ^2(21, t) = in the homogeneous state. Thus 
the first two leading terms i n Eq . ( |12|) can be obtained 
by replacing A 1 + A 2 in Eq. ( |C3| ) with l/2p. 

To obtain the leading contribution to E in Eq. (^) , we 
calculate Ai +A 2 — l/2p, which is expected to be propor- 
tional to dp' 1 /dx. However, the term Ai does not give 
such a contribution. For illustration, it is useful to intro- 
duce a new coordinate x = —x and redefine all quantities 
in terms of the new space variable. Under this transfor- 
mation, p and A\ have even parity [p(x,t) = p(—x,t), 
Ai(x,t) = A\{— x,t)], while the density gradient has 
the odd parity [dp" 1 /dx = —dp~ 1 /dx]. Since A\ and 
dp~ 1 /dx have different parities, A\ should not give a 
correction proportional to dp~ l /dx. 

On the other hand, A 2 gives a correction proportional 
to dp' 1 /dx. One uses the identity 



A 2 {x,t) = ^-[B x {x,t)-B 2 {x,t)] 
dx 



(C5) 



where 



B 1 (x,t) = - J dx'dt'(/)(x - x',t-t') 

x [y r (x',t')(t') - y r (x',t')-i(t')] 2 , 

B 2 (x, t) = 2 dx'dt'(f>{x -x',t- t') 



y r (x',t')(t') + y r (x',t')-i( t ') _ ,' 



(C6) 



In the homogeneous state, B\ 
Thus one obtains 



l/8p 2 and B 2 = l/24p 2 



Here {y n Ay n ) can be approximated by v{Ay n ). Their 
difference is proportional to (dv/dx)(dp~ 1 /dx) (see Ap- 
pendix ^), and thus we ignore (d/dx)[p(y n Ay n ) — 
pv(Ay n )\. One then uses Eq. II to obtain 



(C9) 



By using expansion (12) and the continuity equation ||) 
to convert temporal derivatives into spatial derivatives, 
one obtains 

1 dv 1 d 2 v v ( dp' 1 ^ 
{ Vn) ~~pdx~ + V ^2 + 2 \~dx~ 
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FIG. 1. The velocity profiles for N = 72. The initial con- 
dition in Eq. (Q) is used with A = 1.165 m. (a) t ~ 10 min, 
(b) t w 30 min, (c) t » 1 h, and (d) t « 4 h. The solid 
(dashed) line shows the microscopic (macroscopic) velocity 
profile in each plot. The vertical scale is magnified for clarity. 
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FIG. 2. (a) The time evolution of the space-averaged rel- 
ative deviation of velocity for A = 1.165 m and N — 73. 
(b) fmicro (solid line) vs u maC ro (dashed line) near 115 min 
[marked by the arrow in (a)]. 
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FIG. 3. (a) The density profile for the initial condition in 
Eq. (|3"s|). 8p depends on A and pn- For A — 74.56 m and 
N = 100, Sp ~ 1.5ph- (b) Schematic phase diagrams for the 
microscopic and macroscopic models. 
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FIG. 4. (a) The relative deviations of the macroscopic 
lower (diamonds) and upper (circles) critical densities with 
respect to the microscopic counterparts for the initial con- 
dition [Eq. @] with A = 74.56 m. Note that the relative 
deviations shrink as A increases. (b)The ratio w™ lc /t)™ ac as a 
function of A for the macroscopic model [Eqs. ([]) and (p5h] 
(diamonds) and for the modified macroscopic model [Eqs. (H) 
and (AS)] (circles) takes into account the effects of some lin- 
early irrelevant terms. 



12 



